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Reply: Van den Brink, Khaliullin and Khomskii 

In the preceding Comment Shen points out that 
the on-site Coulomb interaction, that can cause charge 
order in half-doped manganites also destabilizes the 
magnetic CE-phase observed in these systems. 

This is a valid observation, but it is not a priori clear 
whether in the relevant parameter regime the C-phase 
is indeed lower in energy then the CE-phase within our 
model [2). From the exact diagonalization of a 16-site 
ring, Shen finds a groundstate energy for the CE-phase 
as a function of U, which is compared to the energy of 
an infinite ring in the C-phase, in order to find a criti- 
cal coupling Uc — 5t above which the C-phase is stable. 
In order to determine the critical coupling, however, one 
should handle finite size effects in both phases on equal 
footing and compare the CE-phase energy with the en- 
ergy of a 16-site C-phase ring, which is E^q = — 0.628i. 
Using this value for the energy, we can estimate from 
Fig.l in Ref. ||l| that Uc « lOt, which is in the relevant 
parameter regime for the half-doped manganites. 

As Shen states, other interactions which are neglected 
in the model could very well compensate for the small en- 
ergy differences between the C- and CE-phase. We wish 
to point out in particular that coherent and incoherent 
hopping between two sites with opposite spin orientation 
lowers the kinetic energy of the CE-phase appreciably. 

The main aim of our model is to explain the observed 
charge and orbital order, given the CE magnetic struc- 
ture. As to the stability of the phases, indeed not all 
physical relevant effects are taken into account. Espe- 
cially the contribution of lattice distortions, which we ne- 
glected, are important Also, for instance, anharmonic 
effects that favor occupation of elongated orbitals that 
may stabilize the CE-phase, are not incorporated in the 
model. 

Moreover, one should also take into account that the 
hopping between two sites with antiferromagnetically 
coupled core spins is not completely blocked. For a fi- 
nite Hund's rule coupling {Jh) an electron can hop co- 
herently from a high spin state to lower spin states on a 
neighboring site, but even for a infinitely large Jh, there 
exist quantum effects which allow hopping that preserve 
the total high spin state 5-1-1/2, but with different pro- 
jections = S + 1/2 and S — 1/2 at different sublat- 
tices |5] . Such electron hopping processes between C-type 
chains or between CE zigzag chains modify the total en- 
ergy. One can easily estimate the coherent contribution 
to the kinetic energy in second order perturbation the- 
ory. The contribution from hopping between two chains 
in the C-phase is then SEc = — | x and in the CE-phase 
SEcE = 

the CE-phase also the — y 
partially occupied. As typically A w Jh and Jh is in the 
order of 5t, the energy lowering of the CE-phase due to 
this effect would be larger than the energy loss due to the 
on-site Coulomb interactions in the CE-phase. The inco- 


1^. This difference is due to the fact that in 
2 _ ,,2 orbital on the corner site is 


herent contributions that change the spin projections on 
different sublattices are more involved to calculate, but 
they would certainly lower the energy of the CE-phase 
with respect to the C-phase as the inter-chain hopping 
integrals are larger in the former. 

We therefore conclude that the model of Ref. , which 
correctly captures the interplay of spin, charge and or- 
bital degrees of freedom in the half-doped manganites 
and gives a reasonable description of their electronic 
structure, is by itself not sufficient for the precise deter- 
mination of the regions of stabilities of different phases. 
For this several other factors, in particular those men- 
tioned above, should be taken into account. 
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